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Abstract 

In this paper, the quark distribution amphtudes of IS" and 25" heavy quarkonium states are 
studied in terms of Gaussian-type wave functions. The transverse momenta p± integrals of the 
formulae for the decay constant are performed analytically. Then the quark distribution amplitudes 
are obtained. In addition, the ^-moments are also calculated. After fixing the relevant parameters 
appearing in the quark distribution amplitude, the curves of the quark distribution amplitude for 
IS and 2S heavy quarkonium states are plotted. Finally, the numerical results of this approach 
are compared with the other theoretical predictions. 
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I. INTRODUCTION 



Since the discoveries of J /if) and T more than thirty years ago, a great deal of information 
on heavy quarkonium levels and their transitions has been accumulated On one hand, 
the known levels of charmonium and bottomonium include IS", IP, ID, 2S, 2P, 2D, 3S, 
and AS states, with the labels S, P, D corresponding to relative orbital angular momentum 
L = 0, 1, 2 between quark and antiquark. On the other hand, the numerous transitions 
between heavy quarkonium states are classified into strong and radiative decays, which shed 
light on aspects of quantum chromodynamics (QCD) in both the perturbative and the non- 
perturbative regimes (for a recent review see bl). Recently the more precise ip'y Vc^ Vb 
mass measurements have been reported j^, 4, 5, 6|, and the errors of their relevant decay 
widths have decreased Ij. A thorough understanding of their properties, such as their quark 
distribution amplitudes, which are the universal non-perturbative objects, will be of great 
benefit when analyzing the hard exclusive processes with heavy quarkonium production. 

It has been known that heavy quarkonium is relevant for a non-relativistic treatment 

n 

[71]. Although non-relativistic QCD (NRQCD) is a powerful theoretical tool for separating 
the high energy modes from low energy contributions, in most cases the calculation of low 
energy hadronic matrix elements has relied on model-dependent non-perturbative methods. 
Therefore, various methods have been employed in heavy quarkonium physics, such as lattice 
QCD, quark-potential model, etc. j2|]. The light-front quark model offers many insights into 
the internal structures of the bound states. In this study, heavy quarkonium is explored 
within a quark model on the light front. Light-front QCD is a promising analytic method for 
solving the non-perturbative problems of hadron physics j^, and may be the only possible 
method by which the low energy quark model and the high energy parton model can be 
reconciled. For hard processes with a large momentum transfer, light-front QCD reduces to 
perturbative QCD (pQCD) which factorizes the physical quantity into a convolution of the 
hard scattering kernel and the quark distribution amplitudes. 

The basic ingredient in light-front QCD is the relativistic hadron wave function. It 
generalizes the distribution amplitudes by including the transverse momentum distributions, 
and it contains all the information of a hadron from its constituents. The hadronic quantities 
are represented by the overlap of wave functions and can be derived in principle. The 
light-front wave function is manifestly a Lorentz invariant as it is expressed in terms of the 
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internal momentum fraction variables which are independent of the total hadron momentum. 
Moreover, the fully relativistic treatment of quark spins and the center-of-mass motion can 
be carried out using the so-called Melosh rotation [o], |lo|, [ll|. This treatment has been 

ate manyphenomenologically important meson decay constants 
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~v\ . In addition, the covariant light-front 



successfully applied to calcu 
and hadronic form factors 

approach [isl, Q has also been applied to ground-state s-wave mesons, which include 1^5*0 
pseudoscalar mesons r]c,rib and l^Si vector mesons J/ip,T (see Ref. jl8|). 

As mentioned above, the quark distribution amplitude or light-cone wave function 
(LCWF) absorbs the non-perturbative dynamics and is the key ingredient of any hard ex- 



clusive process wit 
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1 hadron production. In the literature, there are many theoretical studies 
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261] of this issue. The main purpose of this study is the calculation 



of the quark distribution amplitudes of pseudoscalar and vector heavy quarkonium states by 
integrating the transverse momenta of momentum distribution amplitudes within the light- 
front approach. As to the so-called ^-moments which parameterize the quark distribution 
amplitudes, they are also calculated analytically. 

The remainder of this paper is organized as follows. Section II comprises brief reviews of 
the light-front framework and the light-front analysis for the decay constants of pseudoscalar 
(P) and vector (V^) mesons; the processes V P7 are given. In Section III, the quark 
distribution amplitudes and the ^-moments are calculated. In Section IV, the numerical 
results and discussions are presented. Finally, the conclusions are given in Section V. 



II. FORMALISM OF COVARIANT LIGHT-FRONT APPROACH 

In heavy quarkonium, the valence quarks have equal masses, mi = 1712 = m, with m the 
mass of heavy quark c or b. Thus, the formulae in this section lead to simplifications for the 
quarkonium system. 

The momentum of a particle is given in terms of the light-front component hj k = 
{k~ , k~^, k±) where fc^ = ± k^ and k_i_ = {k^, k"^), and the light-front vector is written as 
k = {k~^, k±). The longitudinal component k~^ is restricted to positive values, i.e., k~^ > for 
the massive particle. In this way, the physical vacuum of light-front QCD is trivial except 
for the zero longitudinal momentum modes (zero- mode). A meson with total momentum 
P and two constituents, quark and anti-quark whose momenta are pi and p2, respectively 
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will be studied. In order to describe the internal motion of the constituents, it is crucial to 
introduce the intrinsic variables {x,p±) through 

p+ = (l-x)P+, P2± = (1 - x)Px - P±, (1) 

where x is the light-front momentum fraction. The invariant mass Mq of the constituents 
and the relative momentum in z direction can be written as 

The invariant mass Mq of qq is generally different from the mass M of meson which satisfies 
= p"^. This is due to the fact that the meson, quark and anti-quark cannot be on-shell 
simultaneously. The momenta p±_ and pz constitute a momentum vector p = {p±,Pz), which 
represents the relative momenta in the transverse and z directions, respectively. The energy 
of the quark and antiquark ei = 62 = e can be obtained from their relative momenta. 



' m'^ + p^ + pI- (3) 
It is straightforward to find that 



As shown in Ref. [27[, one can pass to the light-front approach by integrating the p~ 
component of the internal momentum in the covariant Feynman momentum loop integrals. 
The Feynman rules for the meson-quark-anti-quark vertices were then needed to calculate 
the amplitudes which related to the decay constant and Ml transition. In the following 
formulations, we follow the notation in [U^. The vertices Tm for the incoming meson M are 
given as 



ifp75 for P, 



Wv 

The Feynman rules which are derived from quantum field theory and Eq. can be used 
to write down the relevant matrix elements. The integration of the p~ component will force 
the quark or anti-quark to be on its mass shell. The specific form of the covariant vertex 
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function for the on-shell (anti) quark can be determined by comparing it to the conventional 
vertex function, which can be written as 

|M(P,2^+i Lj, .],)) = j {(fvi}{d'v2} 2{2Tif5\P -pi- P2) 

X J2 ^^^^(pi,P2, Ai, A2) \qi{pi, Ai)g2(P2, A2)). (6) 
-^1,-^2 

The momentum-space wave function "^^^^ can be expressed as 

^^^^^(pi,P2,Ai,A2) = Rf%{x,p±) (7) 

where (j){x,p±) is the hght-front momentum distribution amphtude for the s-wave meson 
and can be chosen to be normahzable, i.e., it satisfies 



dxd Pj_ 



x,P±)\' = l, (8) 



2(27r)3 

and Rx^x2 constructs a state of definite spin {S, Sz) out of hght-front hehcity (Ai , Ao) eigen- 



states. In practice, it is more convenient to use the covariant form for -Rff^^ 1^' ll^l ]: 



where 



Rx'x.ix, k^) = Ai)rt;(p2, A2), (9) 



r = 75 (pseudoscalar, 5* = 0), (10) 

^ = - ^(^'' + (vector. S = l). (11) 



All details are shown in Appendix A of Ref . [l^ . The function Hpy and the parameter W^ 
are reduced to hpy and wy, respectively, and they are written by 



'V 



/,.^ft.^(M^-M|)^£(l^^^(.,pJ, 

wv = Mo + 2m. (12) 

In principle, 4'{x2,p±) is obtained by solving the light-front QCD bound-state equation 
HlfI"^) = ^1^), which is the familiar Schrodinger equation in ordinary quantum mechanics 
and Hlf is the light-front Hamiltonian. However, except in some simple cases, the full 
solution has remained a challenge. There are several popular phenomenological light-front 
momentum distribution amplitudes that have been employed to describe various hadronic 



structures in the literature. A widely used one is the Gaussian- type chosen here for the IS" 



heavy quarkonium state 



18|: 



3/4 



dpz 



exp 



(13) 



dx \ 2/?2 

In addition, the momentum distribution amplitude for the 25* heavy quarkonium state is 
considered. Eqs. dH]) and ffT^ can be rewritten as 



1, 



3/4 



exp 



2/32 



(14) 
(15) 



where p = \p\- The Fourier transition, or the conjugate coordinate wave function of Eq. ([T5 
is the ground-state solution of the harmonic-oscillator (HO) problem: 

?2 \ 3/4 / n2^2 



-J exp 



Then, for consistency, the relevant excited-state solutions 



IP I 



(3r exp 



3 7ri/4 
1 f(3^^'/' 



/3V 
2 



exp ■ 



(16) 

(17) 
(18) 



Q \7i J V4 2y"V 2 
{YimS are the spherical harmonics), are app lied to the excited meson states. This suggestion 
has been given by the authors of Ref. [28t]. Therefore, the Fourier transforms of Eqs. (flTl) 
and (fT8l) can be rewritten as 

3/4 



(x,px) = 4V2 



TT 



dpz Pm 

dx P 



exp 



pI+pI 



2/52 



0^'^(x,p_l) 



4a/- (- 

''3 



3/4 



dpz 
dx 



pI + pI 

2/52 



4 ) exp 



2/52 



(19) 
(20) 



(p^^^i = =f(p±i ±^P±2)/"\/2, and Pm=o = Pz) and can be treated as the momentum distribu- 
tion amplitudes of the IP and 25* meson states, respectively. In fact, Eq. f|T9l) has been used 
for the p-wave mesons in Ref. 16| . This paper further analyzes the momentum distribution 
amplitudes as shown in Eqs. ( fT3|) and ( l20l) . 

The formulae of the decay constants fpy and Ml transition form factor 1/(0) are needed 



in the latter analysis. However, their derivations have been done in Refs. 
definitions and formulae for them are provided here. 
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18|. The 



A. Decay constant fpy 



The decay constants of mesons fpy are defined by tlie matrix elements for P and V 
mesons 



(0|AJP(P)) = tfpP^, 
{0\V,\ViP,e)) = Mvfve,, 



(21) 
(22) 



where P^ is the momentum of meson and is the polarization vector of V meson. The 
formula for fp is 



47^3 



dx(fp± 



^ ' dx(fp± 



x{l - x){M^ - Ml 

m . , 



m 



<)p[X,p±) 



(23) 



where iY^ = 3 is the color number and m denotes the quark mass. As to the formula for fy, 
we considered the case with the transverse polarization 



and obtain 



V 



e(±) 



Ik 

471-3 



P+ 



e±-P±,0,ei_ 



dx(fp± 



e. = T^(l,±.) 



(24) 



V2K 



dxd'^p±- 



M2 

2 wv 



2 2m 2 



4>v{x,p±). 



(25) 



This expression can be shown to be in agreement with Eq. (2.22) of Ref. [16|]. Since the 
momentum distribution function is even in p^, a quality defined in Eq. ([2]), it follows that 



dxd'p^ '^^^'^'^^ 



1 



X - - ) Mo = 0. 



Therefore, 







dxd'^p± 



vrfi + 2 



(26) 



(27) 



B. Vector current form factor ^(g^) 



For the transition P7, a more general process V — P7* where the final photon is 
off-shell is considered. The V P7* transition is parameterized in term of a vector current 



form factor V^(q'^) by 



(28) 



where is the amphtude of the V P7* process. P (e) is the momentum (polarization 
vector) of the initial vector meson, P' denotes the momentum of the final pseudoscalar 
meson, and the momentum transfer q = P — P'. The formula for the form factor V"(g^) is 



V{q') 



<Pv{x,P±) 
Mo 



(l-x)M^ xM^ 



X 



■2/2 , ip±-q±)' 

m-\ {p±-\ ; 

wv V Q 



(29) 



where p'^ = P± — (1 — x)q±, p'[^ = p± + xq±, Mq = {w? + p'j)/x{l — x), and M, 
{m? + p"l)/x{l — x). The rate for V P7 is 



T{V ^ P7) 



a (M^ - Mf, 



2 \3 



(30) 



III. ANALYSIS OF MOMENTUM DISTRIBUTION AMPLITUDE 

In this section, the momentum distribution amplitudes of pseudoscalar and vector heavy 
quarkonium states are analyzed. The forms of quark distribution amplitude and ^-moments 
can be derived from them. 

A. Quark distribution amplitude ^p{C) 



The quark distribution amplitude of the pseudoscalar heavy quarkonium state can be 
defined as follows [291]: 



Hz)-f'^-f5[z,-z]c{-z)\ 



(31) 



where ^ = 2x — 1 and /i is an energy scale which separates the perturbative and non- 
perturbative regimes. The factor [z, —z] is defined as 



[z, —z] = Pexp[iq J dx^A^^^x)], 



(32) 



which makes the matrix element Eq. (13T!) gauge invariant. The quark distribution amplitude 
$(,^,/x) is normalized as 



y"'rfe$(e,/^)=i, 



(33) 



and it can be expanded 



29| in Gegenbauer polynomials Cn'^{^) as 



$(e,/i) = $a.(0 



n=l 



(34) 



where $as(0 = 3(1 — C^)/4 is the asymptotic quark distribution amphtude and an{lj) the 
Gegenbauer moments which describe to what degree the quark distribution amplitude devi- 
ates from the asymptotic one. Cn'^{^)s have the orthogonality integrals 



Then a„ can be obtained by using the above orthogonality integrals as 

2(2n + 3) 



an(/i) 



cr(0<f(e,^)rfe 



(35) 



(36) 



3(n + l)(n + 2) 

An alternative approach to parameterize quark distribution amplitude is to calculate the 
so-called i^-moments, 



(37) 



It is easy to find the relations between the Gegenbauer moments and the ,^-moments. Here 
we list them up to n = 6 



02 

a4 



ii[i-i4(a+2i(0], 



35 



448 



[5 - 135(^2^ + 495(^^) - 429(^^)]. 



(38) 



All the n-odd moments are vanishing because this paper only focuses on the s-wave momen- 
tum distribution amplitudes, which are £-even functions. 

There are some similar procedures 30|, iMl by which the quark distribution amplitude 
can be obtained from the equal time wave function. Within the framework of this study, 
the decay constant Eq. (1231) can be rewritten as 



dx 



V2K 



and we defined LCWF $p(x,/i) as 



m 



(39) 



<lp(x,^) 



m 



(40) 



where ^p{x,fi) = 2$p(^,/i). If the momentum distribution amphtudes Eqs. ffT^ and (120]) 
are apphed as the non-perturbative inputs, the suppression of the Gaussian function allows 
us to do the p± integrals up to infinity with no loss of accuracy. Thus, the results can be 
obtained as follows: 



, 5/4 

— e'T 


"3 d 


) fF 


4' i_^2_ 



2\ 



5/4 



m 
p 



r 



7 d 



4' 1-^2 



d] r 



d 



where r[r, ?/] is the incomplete Gamma function 



r[r, y] 



t 



'-^e-'dt, 



(41) 

(42) 

(43) 



and d = m?/2(3'^. The incomplete gamma function may be expressed quite elegantly in 
terms of the confluent hypergeometric function 



T[r,y\ = T[r]- r y'' x iFi{r;r + 1; -y), 



where 



iFj{ri,r2, 



r. • ty ly 



{ri)n{r2)n--{ri)n y' 



n=0 



i;nV' 2Jn- 



(44) 



(45) 



and (r)„ = (r+n — l)!/(r — 1)! is the Pochhammer symbol. After fixing the parameters m and 
/?, the behaviors of the quark distribution amplitudes, Eqs. ( HTj) and ( H2l) . are determined 
and compared with those of other theoretical groups in Section IV. It is worth to mention 
that Eqs. fHTj) and (142!) have incorrect asymptotic behavior because the Gaussian-type wave 
functions are primarily only suitable for the low energy region. At energy scales accessible 
at current experiments, a quark distribution amplitude can be far from the asymptotic form. 
The worth of the asymptotic form is that it provides a criterion. In other words, one can use 
the Gegenbauer moments (or ^-moments) to describe to what degree his quark distribution 
amplitude deviates from the asymptotic one. 

In fact, the ^ integrals of $p^^'^(0 can also be analytically performed, and the decay 
constants /p ^''"^ can be expressed as 



15 
P 



25 
P 



5/4 



m e 





'3 
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It is evident that /p^'*'^ is only dependent on the values of m and m/ (3. Eqs. f H6|) and fH7|) 
can be expanded in terms of d as 



IS 



Cl 



~ -c^rf (1 + + + ^rf'-) + c',d''' (1 + 2rf + y rf^ + + ...), (49) 

where all c{ \ positive constants. However, for a typical value of all the series in 
Eqs. f HHj) and fj49l) converge very slowly. In addition, the ^- moments can be expressed 
analytically as 

1 2 



2\15 



15 



P — 



6\15 
P 



lis 



.IS 





"3" 








_4_ 




'3" 








_4_ 




'3" 








_4_ 



iFi 



3 1 
"2' 4' 



2V27r 



7V27r 



ii^i(4i'-^ 



77V27r 



iFi 



, -, —d 

4 ' 4' 



(50) 



and 



7d, 



?4\2S _ a2S \ 



1" 




' 3 






iFi ( 






_"4_ 


V 2' 





(a 



^ - \ 10 
/l5rf 



r 



where 



8V27r 



[-1]^ / 11 3. 



2'4' 



3 /2\^/^ 



308V27r 



2X5/4 



4 '4' 



m 



The derivations of Eqs. (l50l) and (15T1) have used the formula 

6-1 



iFi(a; 6; c) 



LFi(a; 6 - 1; c) - iFi(a - 1; 6 - 1; c)] , 



d] (51) 



(52) 



(53) 



which is easily checked from the definition of the confluent hypergeometric function Eq. 

m. 
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B. Quark distribution amplitude ^viO 



Similar to the case of $p(0; Eq. (125|) can be rewritten and $y(x) defined as 



V2K 



(fp±- 



2m 



Ml 2 — 2 
2 wv 



(54) 



Stt^M fv J ^ + 

However, the integrals in Eq. fl54p cannot be analytically performed . The crux is the 
third term in the square bracket, that is, the term proportional to l/wy- This term may be 
rewritten and expanded as 

2 



1 



1 



Mo + 2m 4m (l + ^^^^^) 



1 

4m 



1 - 



Mo - 2m 
4m 



+ 



Mq - 2m 
4m 



(55) 



because, in the non-relativistic limit. Mo —>■ 2m. Then, we defined the approximate quark 
distribution amplitudes as 

0y(x,p_L) 



V 



V 



d p_L 



■ M2 ■ 

A/m^ + L 2 2 



0y(a;,p±) 



M2 



2 , pi 



Mo - 2m 
4m 



(56) 
(57) 



After applying the momentum distribution function in Eq. f|T3l) and fitting the parameters, 
we found that not only the center values of the parameters (3 of Eqs. ( !56l) and (!57l) were 
both in the error bar of (see Table I), but also that their curves were almost the same. 
In addition, a similar situation also existed in the case of the 25 states. Therefore, we have 
only shown the results from the form ^y{x). By performing the p± integrals, it can be 
obtained that 



$f (0 



\ 5/4 2d 



dV 



d 



2\'l' 



d 



+ (3 + e')r 



7 d 



4' 1_^2 



+ (3 + nr 
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4' l-ei 
d 



(58) 



+ d] idr 



d 



4 ' 1 - 



7 d 



4'1-e' 



.(59) 



The above derivations are based on the formula: 



r[r, y]-ir- l)T[r -l,y]= yiT[r -l,y]-{r- 2)T[r - 2, y]), 



(60) 
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which is easily checked from the definition of incomplete Gamma function, Eq. (H3|) . The 
,^-moments of can also be obtained analytically 



2\1S 
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10 iFi 



3 3 , 
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4V27r 
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4\25 
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V 



where 



425, 



-r 



18 iFi 



32V27r 



7 11 

11 ,F,{--;-;-d 



11 11 

T' T 



15 iFi ( -| -| -rf) + (rf - 6) -d 



3l2 



dV^T [-1] 



425 J 

35^ 



7V27r 
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■4'T'-^, 
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2' 4' 



8\/27r 



35 iFi 



7 11 



Al'i -r 
^ ^ 21 



63 iFi ( -l--^--d] +(rf-46) iFi 



9 3 

'2' -4'-^ 



2 n 



64V27r 



-,(61) 



7 11 \ / 11 11 ^ 

77 iFi ( --; -rf j + 2(d - 46) ( --; -; -d 
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IV. NUMERICAL RESULTS AND DISCUSSIONS 



- 



Try M/f 



(63) 



In this work, the numerical results were calculated for the quark distribution amplitudes 
^{C} and ^-moments for the pseudoscalar and vector heavy quarkonium. First, it was nec- 
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essary to determine the parameters appearing in the momentum distribution functions. In 
total, there are five for charmonium, rric, Pn^, Pj/^p, and five ones for bottomonium: 
mft, Prit, Pri[^, Pr, Pt'- The decay constants of vector heavy quarkonium states were deter- 
mined first. The decay constant /y is related to the electromagnetic decay of vector meson 
V e+e- by [321 



e ' e 



3 M^ 



-cvf^ 



2 



(64) 



where cy is the square of the electric charge of heavy quark. From the data we obtained 
the values 



fj/^ = 416 ± 7 MeV, 
fr = 715 ± 5 MeV, 



/v,' = 298 ± 8 MeV, 
/t' = 497 ± 4 MeV. 



(65) 
(66) 



For the decay constants of rj^, the two decay modes B —>■ Krjc and B —>■ KJ/ ijj are considered 
using the following relation 



r{B ^ Kt], 



T{B^KJ/^) \f. 



Vc 



a 



C 



J/V 



A 



BKric 



X3 

^BK.J/il) 



1 - 



K 



,(67) 



where Cr^^ and Cj/,/, are related to the Wilson coefficients and can be determined to have 



the value = 0.89 ± 0.02 [3 



A 



abc 



1 



M^V ^M^Ml 



1/2 



(68) 



Ml MIJ MlMl_ 
and the form factors f±{q^) are defined by the Lorentz decomposition of the matrix element 



{K{p)\-s^^h\B{p + q)) = U{q'){2p + q), + f^{q^)q,. 



(69) 



By the results of Ref. 16| one calculated the form factors f±{q^) within the covariant 
light-front quark model. Combining the above with the experimental values, we can obtain 



L = 421 ± 38 MeV. 



(70) 



This result leads to the ratio 



1.02 ±0.22, 



(71) 
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which is consistent with the Van Royen-Weisskopf formula 34l |. 









\fj/i'J 







1.04, (72) 

when the approximation "^nd^) — ^j/i/.(0) is used. For the decay constant /,,^, however, a 
similar estimation of Eq. ( 167|) has a large uncertainty because the error of Br{B Kr]'^) was 
greater than 50% [l|. Given the consistency between Eq. flTTj) and fl72l) . the approximation 
\E'^/^(0) ~ \E'^'(0) was used, and the corresponding Van Royen-Weisskopf formula. 



,1) 

is applied and the decay constant /^^ = 300 ± 8MeV is obtained. Furthermore, an ad- 
ditional constraint is needed. The new data Br{J/ip rjc'y) = (1.98 ± 0.31)% |3a] are 
chosen to calculate the form factor V^(0) in Eq. ( l30i) . For the charm sector, there were one 
assumption (Eq. ( |73i) ) and four constraints (the decay constants ( |65l) . (ITOll . and the new 
data Br{J/il) fid)) were used to fix the five parameters appearing in the momentum 
distribution amplitudes. The results are listed in Table 1. As to the decay constants of 

TABLE I: Fixed parameters rric and /?'s for charmonium states (in units of GeV). (5' and (5" are 
the parameters appearing in equations (f56]l and (f57|) . 



parameter 


rric 






l^j/i) 




value 


1.56 


0.820^°:°?! 


0.665l°:°i 


0.613 ±0.006 


0.477 ± 0.008 


parameter 












value 




0.611 ±0.06 


0.613 ± 0.06 


0.474 ± 0.08 


0.477 ± 0.08 



and rj'f^, a similar approximation was assumed 



Jt' J Kit J Mri^ 
and the decay constants fr/^ = 718 ± 5 MeV, /^^ = 499 ± 5 MeV obtained. Finally we used 
[igI i rub = 4.64 GeV and fixed the (3 parameters for the bottom sector in Table II. 

After fixing all parameters, Eqs. fHTl) . (155]) . and (13^ are used to plot the quark 

( "lliS* ( )'2S 

distribution amplitudes (0 and $py (0 for the charm sector, as shown in Fig. 1 and 
Fig. 2, respectively. As is seen, the momentum fraction x of the pseudoscalar charmonium 



15 



TABLE II: Parameters mb and /3's for bottomonium states (in units of GeV). /?' and /?" are the 
parameters appearing in equations ([56|) and 



parameter 


nib 






/3t 


/3x' 


value 


4.64 


1.47 ±0.01 


1.04 ±0.01 


1.30 ±0.01 


0.926 ± 0.005 


parameter 












value 




1.30 ±0.01 


1.30 ±0.01 


0.924 ±0.05 


0.926 ± 0.005 




FIG. 1: Quark distribution amplitudes ^r^dO (solid lines) and ^'j/^iO (dashed lines) of this work. 

had a slightly wider distribution than that of the vector charmonium. Next, the quark 
distribution amplitudes ^-qdO ^i^nd $r?^(0 of this work were compared with those of Ref. 



19|, 



20] in Fig. 3. For the latter, the equation [l9|, 120] 



$(e,/i = 1.2GeV) = c(a,/3,7)(l-e')(«±7e') exp 



(75) 



was used to plot the curves. $^'^(^, /u) and $^'^(^, /i) corresponded to (a = 1, /? = 3.8, 7 = 0) 
jigj ] and {a = 0.027, /3 = 2.49,7 = 1) [3; respectively. In Fig. 3, the major difference 
between the two curves of the IS' states was that the curve of Ref. [l^, on which ^ is peaked 
around zero, was sharper than that of this work. This meant that the momentum fraction 

n 

X in the function used in Ref. [19| was more centered on 1/2 than in the Gaussian- type 
wave function. As to the curves of the 25* state, we found the locations of extreme value by 
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CO 

(N 

e 




FIG. 2: Quark distribution amplitudes ^-q'^iC) (solid lines) and (dashed lines) of this work. 




IG. 3: Quark distribution amplitudes ^r]ciO ^f?c(0 t'^is work (solid lines) and that of Ref. 



19|, 



20( 1 (dashed lines). 



differentiating both Eqs. fH2]) and fl75]) over The results were 



e = o,± 



3 + 4ci' 



e = o,± 



^[37 + /37 - a - (a^ + 2a7 + 6a/57 + 7^ + 6/^7^ + /?V)^^T^^ 



2^7 



respectively. The values are ^ = 0, ±0.463 for this work and ^ = 0, ±0.421 for Ref. [20"]. For 
the bottom sector, the quark distribution amplitudes (0 and (0 a-^e plotted in 
Figs. 4 and 5. These figures show that the x-distribution of the pseudoscalar bottomonium 
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-1.0 



-0.5 



0.0 



0.5 



1.0 



FIG. 4: Quark distribution amplitudes ^nti^,) (solid lines) and <&x(0 (dashed lines) of this work. 



CO 

<N 

e 




FIG. 5: Quark distribution amplitudes ^',j^(0 (solid lines) and $^'(0 (dashed lines) of this work. 



was almost the same as for the vector bottomonium. The reason is that since the differ- 
ences between the pseudoscalar and the vector heavy quarkonium states arise from l/mq 
corrections, the larger the mg, the smaller the difference between them. For comparison, 
the ^-distributions of the charmonium and bottomonium states are plotted together in Fig. 
6. 

In addition, the differences between various quark distribution amplitudes can also be 
shown by ^-moments. We calculated the ^-moments Eqs. ( !50l) . ( ISTI) . ( !6T|) . and ( !62l) . Our 
results and those of other theoretical groups are listed in Tables III and IV. In these tables. 



Refs. [19|, 120, 121I] used the QCD sum rules. The authors of Ref. 



2211 calculated in the frame- 
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-1.0 -0.5 0.0 0.5 1.0 



FIG. 6: Quark distribution amplitudes *^r?b,r)^(C) (solid lines) and ^ric,ri'^{0 (dashed lines) of this 
work. 



work of the BuchmuUer-Tye-potential model and found that their results are in agreement 
with experiments, which included the leptonic widths and hyperfine splittings. The authors 
of 



23| made calculations for the Cornell potential 



are shown in Ref. 



24| . The authors of Ref. 



. The predictions of light-cone sum rules 
25l | not only numerically calculated the 15* 
charmonium states in conventional LFQM, but they also used the variational approach to 
fix the mass of the charm quark. The authors of Ref. 26|] constructed quark distribution 
amplitudes by assuming the very non-relativistic limit, where $(a;) = 6{x — 1/2), and then 
redistributing the parton momenta by relativistic gluon exchange. Depending on the as- 
sumed value of as at the heavy quark scale, they found the ^" moments of rjc in Table III. 
For reasonable values of they obtained qualitatively similar values as in 19| and the 
curves of them are slightly narrower than in our numerical analysis. As to the i^-moments 
for the bottomonium states, they were evaluated and are listed in Table V. 



V. CONCLUSIONS 



This study performed the transverse momenta p± integrals of formulae for the decay con- 
stants fpy of IS' and 25' heavy quarkonium states and then obtained their quark distribution 
amplitudes <l>py(^). In addition, the ^-moments (^^''^'^)p,y were also obtained by integrating 
out C, in $p_y(^). For each heavy quarkonium state, the five parameters uiq, Pp, /3p/, /3y, 
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TABLE III: The ^-moments for the 15 charmonium states. {C^)ric = 



moment 


/C2\ 


{e)j/i' 








/c6\ 


/c6\ 


this work 
flQ 911 


„ -1 -1 „-i-n 010 

n n7n+°-°°'^ 

U.U/'U_QQ07 


0q66+° °°^3 


n n^n^+Q 0052 

0.0307l°:°°2^ 
012+°-°°2 


„ „r.i r+0 0005 

0-02151°:°°°^ 

n ni2+°-°°2 

U.U1^_Q QQ2 


n rvi nn +0 0026 

0.01091°:°°^° 

n032+°-°°°9 

U.UUO^_Q QQQ9 


„ rvn^rrr+O 00023 

0.006571°:°°°^^ 
0033+°-°°°^ 






0.086 






020 




0066 


[23] 




0.084 






0.019 




0.0066 


[2^t 




0.13 






0.040 




0.018 


[25J 


0.084l0:rr 


0.082^°:°°^ 




Q 7+0.001 

"-'•'-'J^' -0.003 


n ni6+°-°°2 

U.U±0_o 002 


nn47+0-0006 


0046+°-°°°^ 

U.UU^D_Q QQJ^Q 


[26j 


0.067 






0.011 




0.004 




TABLE IV: The ^-moments for the 25 charmonium states. (^").^/ = (^")^') 


moment 
















this work 
[2p]t 


Q . 02+0-008 


U.10U_Q QQg 

U--L»_o.07 




0.06001°:°°!^ 


0320+° °°2i 

U.UOZU_Q QQ20 

n nc:i +0.031 
U.UOi_o.031 


0229+°-°°26 
u.u^zy^g QQ23 


00950+° °°°'^^ 
U.UUyOU_g 00074 

p) 7+0.016 
•J-Ui' -0.014 


[22] 




0.16 






0.042 




0.015 


TABLE V: The ^-moments for the 15 and 25 bottomonium states. 


moment 


{eh. 












{eh 


this work 


0.0643i°:r 


0.05981°:°°°! 


0.01031°:°°°} 


0.008941°:°°°}} 


n nn9'?7+0-00005 

U.UUZOr„Q QQQQ5 


noiQ2+° °°°°3 
u.uuiyz_o. 00003 


moment 














ieh' 


this work 


0.08441°:°°!° 


n n79q+o-ooo7 


n ni oQ+0.0003 

U.Ui2»_Qooo3 


0.009841°:°°°}^ 


0.00255l°:°r 


Om 76+0-00005 

u.uui m_o.oooo4 



Pv' which appeared in the momentum distribution amplitude were determined. This study 
first extracted the decay constants fvy from the experimental data of the leptonic decay 
Br(y e^e~), and it then used the Van Royen-Weisskopf formula to obtain the decay 
constants fp^p'. These decay constants were used as constraints to fix the above parameters. 
Then, the curves of the quark distribution amplitudes for the IS* and 25* heavy quarkonium 
states were plotted by the fixed parameters. It was found that, for the IS" charmonium 
state, the momentum fraction x in the function used by Ref. [19|] was more centered on 
1/2 than the one in the Gaussian-type wave function. In addition, the x-distribution of the 
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pseudoscalar bottomonium was almost the same as that of the vector bottomonium. The 
reason for this was that the differences between these heavy quarkonium states, which arise 
from l/mg corrections, become small when itlq is large. Finally, the numerical results of the 
,^-moments were calculated and compared with the experimental data and other theoretical 
predictions. 
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